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ABSTRACT 

Recent studies have shown that strong correlations are observed between the 
low frequencies (1-10 Hz) of quasiperiodic oscillations (QPOs) and the spec- 
tral power law index of several Black Hole (BH) candidate sources, in low hard 
states, steep power-law (soft) states and in transition between these states. The 
observations indicate that the X-ray spectrum of such state (phases) show the 
presence of a power-law component and are sometimes related to simulataneous 
radio emission indicated the probable presence of a jet. Strong QPOs (> 20% 
rms) are present in the power density spectrum in the spectral range where the 
power-law component is dominant ( i.e. 60-90% ). This evidence contradicts 
the dominant long standing interpretation of QPOs as a signature of the ther- 
mal accretion disk. We present the data from the literature and our own data 
to illustrate the dominance of power-law index-QPO frequency correlations. We 
provide a model, that identifies and explains the origin of the QPOs and how 
they are imprinted on the properties of power-law flux component. We argue the 
existence of a bounded compact coronal region which is a natural consequence of 
the adjustment of Keplerian disk flow to the innermost sub-Keplerian boundary 
conditions near the central object and that ultimately leads to the formation of 
a transition layer (TL) between the adjustment radius and the innermost bound- 
ary. The model predicts two phases or states dictated by the photon upscattering 
produced in the TL: (1) hard state, in. which the TL is optically thin and very 
hot ( kT > 50 keV) producing photon upscattering via thermal Componization; 
the photon spectrum index T ~ 1.5 for this state is dictated by gravitational en- 
ergy release and Compton cooling in an optically thin shock near the adjustment 
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radius; (2) a soft state which is optically thick and relatively cold { kT < 5 keV); 
the index for this state, T ~ 2.8 is determined by soft-photon upscattering and 
photon trapping in converging flow into BH. In the TL model for corona the 
QPO frequency v high is related to the gravitational (close to Keplerian) frequency 
i / k at the outer (adjustment) radius and i / low is related to the TL’s normal mode 
(magnetoacoustic) oscillation frequency v MA . The observed correlations between 
index and low and high QPO frequencies are readily explained in terms of this 
model. We also suggest a new method for evaluation of the BH mass using the 
index-frequency correlation. 


1. Introduction 

The short dynamical time scales and coherence of quasi-periodic oscillations (QPOs) 
observed in accreting XRB’s strongly suggest that QPOs can provide valuable information 
on the accretion dynamics in the innermost parts of these systems. In this respect the 
discovery of low 20-50 Hz QPOs in luminous neutron star (NS) binaries by van der Klis 
et al. (1985), kilohertz QPOs in NS’s by Strohmayer et al. (1996) and hectohertz QPO’s 
in BH’s by Morgan, Remillard k Greiner (1997) opened a new era in the study of the 
dynamics near compact objects. Following these discoveries Psaltis, Belloni & van der Klis 
(1999), hereafter PBK, demonstrated that these NS and BH low and high frequencies follow a 
remarkably tight correlation. These features are the “horizontal branch oscillations” (HBO) 
along with “low frequency noise Lorentzian” , viow and the lower kHz QPO Vhigh f° r NS and 
BH respectively. Belloni, Psaltis k van der Klis (2002), hereafter BPK, have updated 
PBK’s correlation adding data from Nowak (2000), Boirin et al. (2000), Homan et al. 
(2001), Di Salvo et al. (2001) and Nowak et al. (2002). PBK suggest that the low and high 
frequencies correlate in a way that seems to depend only weakly on the properties of the 
sources, such the mass, magnetic field, or possibly the presence of a hard surface in compact 
object. Mauche (2002) has reported low and high frequency QPO’s in the dwarf nova SS 
Cygni and VW Hyi (see Wouldt k Warner 2002 for for the VW Hyi observations) and has 
called attention to the fact that these QPO’s extend the correlation of PBK downward in 
frequency by more than two orders of magnitude This frequency correlation observed over 
such a broad range of sources indicates that a common phenomenon is responsible for QPOs 
in compact objects and raises the question whether one can find a observational relations 
which would distinguish BH candidates from other compact sources. 

The comparative study of spectral and QPO features in BH candidate sources has 
yielded important information which can help answer this question. As far as spectral 
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properties are concerned, there is increasing observational evidence that the large complex 
of spectral “states” originally devised by Belloni (2000) to classify the highly variable X- 
ray spectra of BH galactic sources like GRS 1915+105, can be reduced to a few simple 
configurations [see e.g. McClintock and Remillard (2003) and Belloni (2003)]: (1) a low 
hard power law dominated state where the photon spectral index is near 1.5; (2) a steep 
power law state in which the index is about 2.7; (3) and a thermally dominated state. The 
observed variability is explained by transitions between these three canonical configurations. 

Low frequency (1-10 Hz) QPO’s are generally observed to present in states (1) and (2), 
during transitions between states (1) and (2) but never in (3). Both high and low frequency 
QPO’s have been thought to be associated with the accretion disk, though their specific 
origin has not widely understood (although see Titarchuk & Wood 2002, hereafter TW02, 
for explanation of their origin). Attempts to find consistent correlations with observable 
disk parameters have lead to mixed results. In particular, there are so many exceptions 
to the observed correlations of frequency to disk spectral values sometimes observed that 
no general statement about the relation of the QPO with any spectral parameter solely 
associated with the disk, such as thermal flux, temperature, inner radius, etc. can be made. 
For example, exceptions to one of the most widely quoted correlations, i.e. between disk 
flux and frequency, have been presented for the well studied micro-quasar GRS 1915+105 
by Fiorito, Markwardt, & Swank (2003) and for XTE 1550-564 by Remillard et. al. (2002). 

On the other hand, a number of recent studies have revealed consistent and robust 
correlations between the photon index of the power-law component of the X-ray spectrum 
and low QPO frequency. Such strong correlations are observed the well-studied BH binary, 
GRS 1915+105 in almost all situations in which the QPO is observable and in several other 
BH’s over a wide range of observations and states. 

These studies include: 

a. A comprehensive study of correlations between photon index and QPO frequency in 
BH sources by Vignarca et al. (2003); these results show strong correlations between index 
and QPO frequency as well as saturation of the index between at the values 1.5 ± 0.1 and 
2.7 ± 0.2 for several BH sources and for the same source in different states. 

b. Recent studies by Kalemci (2002) for a number of galactic BH x-ray transients 
observed during outburst decay which confirm the persistence of the correlation between 
frequency and index during state transitions. 

c. Our own studies of GRS 1915+105 of x states (Fiorito, Markwardt, & Swank 2003), 
in which strong correlation of QPO frequency with power-law index have been observed; 
such states have been traditionally noted as low brightness, power law dominated states (see 
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e.g. Belloni et al. 1999) where the disk is at least tenuous if not undetectable altogether. 

d. Studies tracing the QPO frequency and index over time between state transitions 
shows that the frequency tracks the index during state evolution (Homan, et. al. 2001) and 
Sobczak et al. (1999, 2000) from a value of about 1.5 ± 0.1 to about 2.7 ± 0.2, as the QPO 
frequency varies from about 1 to 10 Hz. 

The persistence of the correlation of index with QPO frequency and its tracking with 
respect to time suggests that the underlying physical process or condition which gives rise to 
the low frequency QPO are tied to the corona, which gives rise to the power-law component, 
and furthermore, that this physical process or condition varies in a well defined manner as the 
source progresses from one state to another. Moreover, the fact that the same characteristics 
of the correlation are seen in so many galactic X-ray binary BH sources, which vary widely in 
both luminosity (presumably with mass accretion rate) and state, suggests that the physical 
processes controlling the index and the low frequency QPOs is a characteristic property of 
these sources, and by virtue of the observed correlations of PBK, may be a universal property 
of all accreting compact systems. 

The above data have motivated us to develop a detailed model of the physics of the 
corona surrounding a BH s which directly predicts the behavior of the spectral index with 
fundamental properties of the corona. The model we have developed incorporates funda- 
mental principles of fluid mechanics, radiative transfer theory and oscillatory processes. It 
identifies the origin of the QPO as a fundamental property of a compact coronal region 
near the BH and shows how the photon index of this corona changes as a function of mass 
accretion rate. 

It has been already shown in a number of papers (see below) that the absence of the 
firm surface in the BHs leads to the development of the very strong converging flow when 
the mass accretion rate is higher than Eddington rate Msdd — Lsdd/c 2 - The main observa- 
tional features of the converging flow should be seen in high/soft phase, while the thermal 
Comptonization spectrum of the soft (presumably disk) photons should be seen in the hard 
phase. 

In series of papers (Chakrabarti & Titarchuk 1995, hereafter CT95; Titarchuk, Mas- 
tichiadis & Kylafis 1996, 1997; Ebisawa, Titarchuk & Chakrabarti 1996; Titarchuk & Zannias 
1998; Shrader & Titarchuk 1998, 1999; Borozdin et al. 1999; Laurent & Titarchuk 1999, 
2001, hereafter LT99 and LT01 respectively; Titarchuk & Shrader 2002; Turolla, Zane & 
Titarchuk 2002) the authors argue that the drain properties of black horizon are necessarily 
related to the bulk inflow in BH sources and that the spectral and timing features of this bulk 
inflow are really detected in the X-ray observations of BHs. The signatures of the inflow are: 
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i. an extended steep power law where the photon index saturates to the value 2. 7 ±0.2 as the 
mass accretion rate increases - the precise value of the index is a function of the temperature 
of the flow but its relatively high value is a result of inefficient photon upscattering in the 
converging flow due to photon trapping. Moreover, LT99 demonstrated using Monte Carlo 
simulations that in the low hard state when the plasma temperature is of order 50 keV the 
bulk inflow spectrum is practically identical to thermal Comptonization spectrum (there is 
no any noticeable effect of the bulk Comptonization in the spectrum) and the spectral index 
varies in narrow range from 2 to 1.7 with mass accretion rate (in Eddington units) increases 
from 2 to 7. The similar small variation of the photon index around 1.7 in the low hard state 
was pointed out in earlier work by CT95. In fact, the effect of the bulk Comptonization 
compared to thermal one is getting stronger when the plasma temperature drops below 10 
keV. 


ii. the QPO high frequency (100-300 Hz), which is inversely proportional to BH mass; 
iii these two features are closely related but the variability of the soft spectral component is 
not correlated at all with the variability of the power law component. 

Furthermore, Titarchuk, Osherovich fc Kuznetsov (1999), hereafter TOK, presented 
observational evidence that the variability of the hard X-ray radiation in NS’s and BH’s 
occurs in a bounded configuration. Ford & van der Klis (1998) found that the break frequency 
in the power density spectrum (PDS) is correlated with QPO frequency for some particular 
NS sources (4U 1728-32). Wijnands & van der Klis (1999) later found a similar correlation 
in BH sources. Titarchuk & Osherovich (1999), hereafter T099, and TOK explained t his 
correlation in NSs and BHs respectively. They argue the observed PDS is the power spectrum 
of an exponential shot, exp {—t/td) which is the response of the diffusion propagation to 
any perturbation in a bounded medium. The inverse of the characteristic diffusion time 
Ub = 1/td and the normal mode QPO frequency v 0 has a well defined relation to the size 
of the bounded configuration, because u b oc V/L and u 0 oc V/L 2 where V is the specific 
pertrubation (magneto-sonic) velocity and L is a charactestic size of the configuration (see 
T099). 

Titarchuk, Lapidus k, Muslimov (1998), hereafter TLM98 proposed that this bounded 
configuration (cavity) surrounding compact objects is the transition layer (TL) that is formed 
as a result of dynamical adjustments of a Keplerian disk to the innermost sub-Keplerian 
boundary conditions. They argued that this type of adjustment is a generic feature of the 
Keplerian flow in the presence of the sub-Keplerian boundary conditions near the central 
object and that it does not necessarily require the presence or absence of a hard surface. 
TLM98 concluded that an isothermal sub-Keplerian transition layer between the NS surface 
and its last Keplerian orbit forms as a result of this adjustment. The TL model is general and 
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is applicable to both NS and black hole systems. The primary problem in both NS and BH 
systems is understanding how the flow changes from pure Keplerian to the sub-Keplerian as 
the radius decreases to small values. TLM98 suggested that the discontinuities and abrupt 
transitions in their solution result from derivatives of quantities such as angular velocities 
(weak shocks). They were first to put forth the possibility of the transition layer formation 
to explain most observed QPOs in bright low mass X-ray binaries (LMXBs). In Figure 1 we 
illustrate the main idea of the transition layer concept for NS’s ans BH’s. 

Subsequently, T099 identified the transition layer (TL) model’s corresponding viscous 
frequency vy and diffusion or break frequency V\, for 4U 1728-34. The predicted values for 
v b related to diffusion in the transition layer. Then, TOK confirmed, with observations of a 
number of QPO sources (Wijnands et al. 1999), the relationship v b a z/^ 6 , which is valid for 
both NS and BH candidate sources. 

In this Paper, we explain the general correlation between photon index and low frequency 
in the frameworks of the transition layer model and we give more arguments for the nature 
of the spectral phases (states) and phase transition observed in BHCs. The main features of 
the TL model are given in §2.1. The formulation of the problem of low frequency oscillations 
and the relationship with MA oscillations in the transition layer are described in §2.2. The 
derivation of the low frequency-index correlation is present in §2.3. We give the details of 
modeling of the spectral phase transition in §2.4. In §3 we present and analyze specific QPO 
and spectral data in terms of the TL model. We discuss the signatures and the methods of 
the identification of BH and NS sources using timing and spectral characteristics in §4. Our 
summary and conclusions also follow in §4. 


2. Transition Layer Theory 
2.1. The main features of the transition layer model 

TLM98 define the transition layer as a region confined between the the inner sub- 
Keplerian disk boundary and the first Keplerian orbit (for the TL geometry, see Fig. 1 and 
Fig. 1 in TLM98 and TOK). The radial motion in the disk is controlled by friction and 
the angular momentum exchange between adjacent layers, resulting in the loss of the initial 
angular momentum by accreting matter. The corresponding radial transport of the angular 
momentum in a disk is described by the following equation (see e.g. Shakura & Sunyaev 
1973, hereafter SS73): 

M^R 2 ) = 2 n—(W rr R 2 ), 


( 1 ) 


-7- 


where M is the accretion rate in the disk and W Ttp is the component of a viscous stress tensor 
that is related to the gradient of the rotational frequency u> = 27 tv, namely, 

W rr = -2 nHR~, (2) 

where H is a half-thickness of a disk and 77 is turbulent viscosity. The nondimensional 
parameter that is essential for equation (??) is a parameter 7 (often called the Reynolds 
number) which is the inverse of the SS73 a parameter for the accretion flow, 

7 = M/farqH = RV r /D, (3) 

where V r is a characteristic velocity, and D is the diffusion coefficient. D can be defined 
as D — V t l t /3 using the turbulent velocity and the related turbulent scale, respectively or 
as D = Dm = c 2 /a for the magnetic case where o is the conductivity (e.g. see details 
of the D-definition in Lang 1998). We present u>(R) in terms of dimensionless variables, 
namely , angular velocity 9 = u/u 0 , radius r = R/Rq (Rq — x 0 Rs, and Rs = 2 GM/<? is the 
Schwarzchild radius) and mass m = M/M & . We express the Keplerian angular velocity as 

9 k = 6/(a K r 3/2 ), (4) 

where ok = m(x 0 /3)(i'o/363 Hz); the rotational frequency u 0 has a particular value for each 
star. Equations 9 — 1 at r = rj n (the inner disk radius), 9 — 9k at r = r mt (the radius 
where the transition layer adjusts to the Keplerian motion), and d9/dr = d9 K /dr at r — r^t 
were assumed by TLM98 to be the boundary conditions. Thus, the profile 9{r) and the 
outer radius of the transition layer are uniquely determined by the boundary conditions and 
the angular momentum equation (1-3) for a given value of 7 parameter. The solution of 
equations (1-3) satisfying the above boundary conditions is equation (10) in TLM98 and the 
equation (see TLM98, Eq. 11) 

3<W2 = At rXk + 2(1 - A)r- 2 , (5) 

determines r^, = R^j R q as a function of 7 parameter, where D[ - (O^i — r «3)/( r «3-0- 
The adjustment of the Keplerian disk to the sub-Keplerian inner boundary creates conditions 
favorable for the formation of a hot plasma outflow at the outer boundary of the transition 
layer (TLM98), because the Keplerian motion (if it is followed by sub-Keplerian motion) 
must pass through the super-Keplerian centrifugal barrier region. For the observed range 
from 1 to 100 Hz and from 0.1 to 10 Hz for of the high and low QPO frequencies respectively 
(see Fig. 2) one can find using the TLM98 solution for m=10 that 7 varies from 2 to 20 (a 
variation which is related to the mass accretion rate M (see Eq. 3). 

It is worth noting that in the literature that a — is treated as a constant and a free 
parameter independent of the mass accretion rate (or the spectral state), whereas we find 
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that the a varies from 0.5 to 0.05 with M, i.e. a is a real strong function of M. This a 
variation has many consequences for the accretion disk model as well as in the interpretation 
of data (see discussion in Titarchuk & Shrader 2002). 


2.2. Magnetoacoustic oscillations in the disk transition layer 

The TL model identifies the low frequency QPO as that associated with the viscous 
magnetoacoustic (MA) oscillation of the bounded TL; this mode is common to both BH’s and 
NS’s. The correlation of the MA frequency with the gravitational (Keplerian) frequency v K = 
2ttlo k were derived by T099, Titarchuk, Bradshaw & Wood (2001) for NS and then it was 
generalized by TW02, for BH and WD. The MA frequency is derived as the eigenfrequency 
of the boundary-value problem resulting from a MHD treatment of the interaction of the 
disk with the magnetic field. The problem is solved for two limiting boundary conditions 
which encompass realistic possibilities. The solution yields a velocity identified as a mixture 
of the sound speed and the Alfven velocity. The TBW treatment does not specify how the 
eigenfrequency is excited or damped. However, it makes clear that the QPO is a readily 
stimulated resonant frequency [see Titarchuk, Cui Wood (2002) and Titarchuk (2002) for 
details of excitation of the eigenfrequencies] . 

A linear relation was derived between v MA and u K (see Fig. 2 and TW02) 

vma — Cmav k- (6) 

where C M a = 2 1/2 47 t [(1 + //?)/( 1 + /5)] 1/2 (F/r out ) is a proportionality coefficient which is 
“universal” to the extent that f3 and HjR 0U t remain about the same from one source to the 
next. Here (3 = P 9 /Pb is the ratio of the gas pressure to the magnetic pressure and the 
coefficient / = 1 /32, 1 /327T 2 for stiff and free boundary conditions respectively and H is the 
half-width of the Keplerian part of the disk. 

Identification of v high , v low , with u K , vma respectively, leads TW02 to determination of 
H/Rout = 1-5 x 10 -2 and /? = 10 using the observed value of Cma = 0.08. 

In Figure 3 we present the inferred v low (u iow = 0.08 v high ) as a function of y. We derive 
this relation using equations (4) and (5) with the assumptions that an = 771(^0/363 Hz) = 1 
and m = M/M 0 = 10. 

It is evident that the inferred Vhigh, 11 low are inversely proportional to m because 

v high = VUL OC r/T 1 (7) 


and 


viow ~ V/L, 


( 8 ) 
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where L is the TL characteristic size is proportional to m. 

We calculate the size of the TL, r mt = Rout/Rs as a function of using relations (6), 
and Vhxgh — ^k- Because the radius of interest, which in our case is the outer TL boundary 
radius Rout, is related to u K one can relate R^t to v iow using formula (6). We present this 
function in Figure 4. 

The spectral index— frequency correlation can be derived if we find a spectral index— TL 
optical depth correlation. Because 7 is proportional to the mass accretion rate M and the 
optical depth tq is also proportional to M , we are lead to a correlation of ui ow and r 0 and 
ultimately to the index— frequency correlation. 


2.3. Spectral index-optical depth relation. Corona model 

It is an worthwhile and interesting question to ask if there is a common mechanism in the 
energy release (heating) in coronae of the compact objects (neutron stars and black holes). 
It is known from X-ray observations of BHs and NSs the emission from the hot coronal 
plasma (source of the hard radiation) is a significant fraction of the bolometric radiation. 
Particularly in the low/hard state the hard radiation (Comptonization) is dominant 5 . In 
hard state spectra in BHs and NS show are both well fitted by power law spectra which is 
a signature of Compton upscattering of the soft disk photons. The BH hard state spectra 
are characterized by spectral indices a = 0.6 ± 0.2 (or photon indices T = 1.6 ± 0.2) and the 
plasma temperature about 50-100 keV independently of the luminosity level. In this state 
thermal Comptonization is dominant. 

In this paper we argue that the heating of corona (Compact cloud) can be determined 
by the gravitational energy release there. 

TLM98 offered the TL model adjustment of the Keplerian accretion flow (disk) to the 
sub-Keplerian innner boundary condition which provides a compact corona for Comptoniza- 
tion upscattering of the illuminating accretion disk photons. The essential sub-Keplerian 
rotation of NS in LMXBs is a well established observational fact due to the direct mea- 
surements of NS spins (see a review by van der Klis 2000 for details). The direct BH spin 
measurements are still unavailable but if we use analogy with NS rotation we can assume 
that BH may also undergo sub-Keplerian rotation. 


5 As the energy release in the solar corona 1.4 x 10 28 erg s, 1 is only small fraction of the solar irradiance 
3.8 x 10 33 erg s' 1 . 
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It is likely that this adjustment to sub-Keplerian flow is not smooth and that near the 
adjustment radius the strong or weak shocks can be formed. For a given viscosity 77 the 
position of the adjustment is determined by the disk mass accretion only (see Eqs. 3 and 5). 
The adjustment radius (shock) region can be treated as a potential wall at which the accreting 
matter releases its gravitational energy in a geometrically thin target. As a consequence of 
this the shock plasma temperature is much higher than that in the surrounding regions of 
the Keplerian disk. The shock region gets puffed up oscillating with frequencies close to u K 
(Titarchuk 2003). Additional oscillations can appear (at least in NS case) because of the 
rotational configuration above the disk which is a magnetosphere in NS case (see details in 
TOK). The shock formation leads to the formation of the hot inner bounded region. The 
hard photons in the shocked region illuminate the surrounding material evaporating some 
part of the disk. It is obvious that the evaporated fraction depends on the mass accretion 
rate in the disk (see Chakrabarti & Titarchuk 1995, hereafter CT95). Some small fraction 
of the hard photons (at maximum 6 %) can be reflected by the cold parts of the disk (CT95 
and see also Basko, Sunyaev & Titarchuk 1974). 

Let us assume the column density of this region y 0 is a few grams or that the Thomson 
optical thickness r 0 is a few (see SS73 for the thickness estimate of the innermost part of 
the disk). The rate of energy release at this region Q cor is a few percent of the Eddington 
luminosity since the adjustment radius is located within 4-100 Schwarzschild radii (see Fig. 
4). The heating of gas due to the gravitational energy release should be balanced by the 
photon emission. For high plasma temperature, Comptonization and free-free emission is 
the main cooling channel, and the heating of electrons is presumably due to the Coulomb 
collisions with protons (see Fig. 1 for the picture of disk corona). 

Under such physical conditions the energy balance can be written as [see also Zeldovich 
& Shakura (1969), hereafter ZS69 and TLM98] 

Ocor/To = <W(T)r.//(r.) + c„TWp, (9) 

where r is the optical depth within the shock, e(r) is a distribution function for the radiative 
energy distribution, f(T e ) = 1 + 2.5 (kT e /mc 2 ) is the relativistic correction factor, T e is the 
electron temperature in K, C comp = 20.2 cm s _1 K _1 and Cff = 2.6 x 10 20 erg K -1 / 2 cm s -1 
are a dimensional constants. In this formula we neglect the gas heating due to recoil effect. 

The distribution e(r) is obtained from the solution of the diffusion equation (ZS69), 

d 2 e/dr 2 = -3(Q tot /c)/r 0 (10) 

subject to two boundary conditions. Here we assume that the region of the gravitational 
energy release (corona) is a spherical shell surrounding the central object and the total 
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flux in the corona is a sum Q tot = Qcor + Qd of the gravitational energy release Q C0T and 
the illumination flux from outside of the corona (disk) Q d . The inner and outer boundary 
conditions for a BH corona are that there are no scattered radiation from outside of the 
corona. Then in the Eddington approximation (see Sobolev 1975) the boundary conditions 
can be written as follows : 

(de/dr - 3s/2)| T=0 = [del dr + 3e/2)| T=To = 0, (11) 

where r = 0 and r = tq are at the inner and outer boundaries respectively. 

For a NS corona one should assume that inside of the inner corona boundary the energy 
density e =constant, namely that the radiation flux emitted towards the central object at 
some point of the inner boundary returns back (reflected by the NS surface), namely this 
condition is equivalent to the reflection condition: 

de/dr\ T=0 = 0. (12) 

Below we present a general formula of the coronal temperature which combines the BH and 
NS cases. 

The solution of equations (10-11) provides us with the distribution function for the 
energy density in the BH corona 

£ ( T ) = (Qtot/c){ 1 + (3/2 )t 0 [t/t 0 - (t/t 0 ) 2 ]}. (13) 

We neglect the dependence of Q to t on r to derive this formula. Furthermore, we take a 
representative value of < e > for the maximum value in the corona which is very close to 
mean value of e mean i.e. 

< £ Bmax = (Qtot/c)[ 1 + (3/8) To]. (14) 

in order to approximate the coronal temperature. In Figure 5 we present the results of calcu- 
lations of the temperatures kT e (keV) and the energy spectral indices of the Comptonization 
spectrum a (photon index T = a + 1) asa function of the optical depth of the shell To. 
The different curves of T e (r 0 ) are related to the different ratios of Qd/Qcor • The calculations 
have been made for values of p = 10 -6 g cm -3 and — 10 22 erg cm -2 s _1 which are 
characteristic values of the density and luminosity for the standard disk model (SS73). 

In fact, the calculation results weakly depend on p and Q C0T if 
(p/10 -6 g cm -3 )/ (Q cor /10 22 erg cm -2 s -1 )^l 


is of order a few or less. 
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In this case the equation for the temperature is simplified 

6.73 x 10 -2 to(1 + 3r 0 /8)(T e /10 8 K)( 1 + Qd/Qcor)/ f(T e ) = 1 (15) 

which has the solution 

r e /10 8 K = 1.5 x [r 0 (l + 3ro/8)(l + Q d /Q cor ) - 0.62]- 1 . (16) 

It is easy to show that the case for which Qd/Qcor = 1 is identical to the NS case (the inner 
reflective boundary) when the disk illumination is neglected. One can see from Eq. (16) 
that T e is a function of r 0 if Qd/Qcor < 1- It is also worth noting that T e r 0 ( 1 + 3r 0 /8)//(T e ) 
is insensitive to the total luminosity Q to t if if Qd/Qcor <C 1 (Eq. 15). 

The spectral index aasa function of the product of T e r 0 are insensitive to Q tot too for 
1 because the the Comptonization parameter 

y = 4kT e /m e c 2 Max(r 0 ,ro) 


and 

a = -3/2 + V(9/4 + 4 Jy) 

(see Rybicki & Lightman 1979 and Sunyaev & Titarchuk 1980). This inferred property of 
the Comptonization spectra reproduces the observed independence of a on the bolometric 
luminosity (see e.g. Tanaka 1995) when a varies within the range 0.6 ± 0.2 (r = 1.6 ± 0.2). 
This defines the so called the low/hard state in BHs. 

Below we show that the observed plateau of the index-frequency correlation at low values 
of v L < 1 Hz is the result of this behaviour because the lower frequency values are related 
to the low mass accretion rates m = M/Msdd < 1 and r 0 is of order one when rh < 1. 

We calculate the exact values of the spectral indices a for a given T e and r 0 using 
the relativistic formulas developed by Titarchuk & Lybarskij (1995) [see formulas (17), (for 
plane geometry, note their r 0 is our tq/2) and (24) there]. In Figure 5 we present the spectral 
indices as a function of r 0 for BH and NS cases when the disk illumination is negligible small 
with respect to the coronal energy release (a condition which can be a model for the low/hard 
state in these systems). The plasma temperature values of 20-150 keV for r 0 = 1 — 3 and 
a = 0.6 ± 0.2 (for ratio = Qd/Qcor — 0) are typical values of these quantities for low/hard 
state in BHs. 

When Qd is comparable with Q cor the cooling becomes more efficient due to Comp- 
tonization and free-free processes and therefore T e unavoidably decreases [see T e — values for 
Qd/Qcor = (1 -4) in Fig. 5]. 
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This illumination effect and mass accretion rate increase may explain the hard-soft 
transition when the Compton temperature drops substantially with increase of the soft (disk) 
photon flux (for illustration, see Figure 1 ). In Figure 5 we show that temperature drops 
from 50-60 keV for ratio — 0 to 7-10 keV for ratio = 4. The temperature may drop also due 
to an increase in the optical depth (presumably because m increases) even when ratio = 0. 


2.4. Modeling of spectral phase transition from the hard state to the soft state 

Laurent and Titarchuk (1999), hereafter LT99, studied the Comptonization of the soft 
radiation in the converging inflow (Cl) into a black hole using Monte Carlo simulations. 
The full relativistic treatment has been implemented to reproduce the spectra. An accreting 
black hole is, by definition, characterized by its drain effect , namely, matter falls into 
a black hole much the same way as water disappears down a drain: matter goes in and 
nothing comes out. The effect of the drain within the TL corona, i.e. the converging inflow, 
bn the X-ray spectrum provides a unique observational signature of black holes. LT99 show 
that spectrum of the soft state of BHs can be described as the sum of a thermal (disk) 
component and the convolution of some fraction of this component with the Cl upscattering 
spread (Green’s) function. The latter boosted photon component is seen as an extended 
power law at energies much higher than the characteristic energy of the soft photons. LT99 
demonstrate the stability of the power law index (the photon index, T = 2.8 ± 0.1) over a 
wide range of the plasma temperature 0-10 keV, and mass accretion rates (higher than 2 in 
Eddington units) due to upscattering and photon trapping in the converging inflow. The 
spectrum is practically the same as that produced by standard thermal Comptonization (see 
e.g. Titarchuk 1994, Hua & Titarchuk 1995) when the Cl plasma temperature is of order 50 
keV (the typical ones for the BH hard state). In this case one can see the effect of the bulk 
motion only at high energies, where there is an excess in the spectrum with respect to the 
thermal Comptonization spectrum. 

LT99 also demonstrate that the change of the spectral shapes from the soft state to the 
hard state is clearly related to the temperature and optical depth of the bulk inflow. We 
combine our results of calculations of index and Compton cloud temperatures for thermal 
Comptonization case (see §2.3 and Fig.5) and LT99’s results of the spectral calculations (see 
Table 2 in LT99 for the spectral index values) to evaluate the power law index as a function 
of the optical depth of the Compton cloud r 0 . We use the values of the index and of the 
temperature for Q = 0 (see Fig. 5) when To varies from 1 to 3 to decribe the index behavior 
and the plasma temperature variation in the hard state. Figure 6 presents the photon index 
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r as a function of r 0 . This function r(r 0 ) can be fitted analytically as follows: 

r(r 0 ) = c 0 4- Ci To + c 2 Tq + c 3 Tq for 1 < r 0 < 5.17, 

r(r 0 ) = b 0 - 5ie- To/T * for r 0 > 5.17 (17) 

where c 0 = 0.61, c x = 0.57, c 2 = 1.485 x 10" 2 , c 3 = -8.5 x 10~ 3 and b 0 = 2.81, b x = 1.0 x 10 3 
and r* = 0.5. 

In order to derive the index dependence V on ui ow one should relate the parameter 7 to 
r 0 and then combine the two relationships T vs t 0 and vs 7 (see Figures 3-6 ). In the 
first approximation t 0 should be proportional to 7 because both of them are proportional 
M. To fit the data we assume that t 0 = A'y 8 . We have three free parameters, black mass 
m = M/M©, 5 and A to fit the model to the data if we assume that the converging inflow 
(CI) is essentially cold (the CI temperature kT c i < 5 keV). 


3. Data and Analysis 

We analyze data points for photon index vs low frequency correlation obtained from 
Trudolyubov et al. (1999) and Vignarca et al. (2003) for GRS 19154-105 and Sobczak 
et al. (1999, 2000) for XTE J1550-564. Also we make use data for GRS 19154-105 by 
Fiorito, Markwardt & Swank (2003) hereafter FMS, which is a detailed study of QPO and 
spectral features of the steady low hard states of GRS 19154-105 for representative 1996 
and 1997 observations. FMS computed the power spectra of combined Binned and Event 
Mode data for PCA channels 0-255 and sampled at 0.25-64 Hz for 4-second intervals. The 
spectra were averaged over the entire observation duration. The PDS spectra were fit with 
a combination of Lorentzian lines and exhibited a break frequency, a central QPO frequency 
and in some case the presence of the first harmonic. Background subtracted energy of each 
OBSID were obtained and time average over the entire observation duration using Standard 
2 PCA and HEXTE data where available (both clusters). The XSPEC model used to fit 
the energy spectra consisted of a multi-temperature thermal disk component (DISKBB), a 
simple power law and a small Iron line component. For the 1996 data (OBSID10408), it was 
found that the disk component could be neglected without affecting the fit. Acceptable y 2 
were obtained in all cases. 


3.1. Correlations between QPO Frequency and Spectral Parameters 

The first part of the FMS study was to first to look for correlations of different spectral 
components, presumably thermal disk components, with QPO frequency to try to determine 
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the reason for the non-unique total flux correlations previously observed. In the case of 
OBSID 20402 (1997) we did observe the expected positive correlation of QPO frequency 
with disk flux in agreement with previous analyses of Muno et al. (1999) and Markwardt 
et al. (1999), which were done in other states. However, we observed an anti-correlation 
of frequency with disk flux for OBSID 10408 (1996); and no correlation of frequency with 
disk flux for OBSID 10258 (1996). Furthermore, no correlation of QPO frequency with disk 
temperature is seen for any of the OBSIDS we studied. These results are at variance with 
the often-quoted result that disk parameters are closely associated with QPO frequency. 
Remillard, et. al. (2002) show a comparison of low QPO frequency and the apparent disk 
flux for the source XTE J1550-564 during its 1998-1999 outburst. The QPOs are reported 
with differences in phase lags. Again, it is clear that a frequency-disk flux correlation is 
observed for this source is only observed at lower (v < 5 Hz) QPO frequencies and fluxes 
than usually observed for GRS 1915+105. 

However, robust correlation between index and the central frequency of low frequency 
QPOs are also seen in a wide variety of BH sources radiating in different states and during 
transition between states. The QPO frequencies in such studies are observed in the PDS 
obtained by integrating over time intervals ranging from the complete time of the observation 
to segmented time bins where the luminosity remains constant, [see e.g. Markwardt, et. al 
(1999)]. Strong frequency-index correlations are observed by Kalemci (2002) in outburst 
decay (presumably state transitions) of a number of sources and by Vignarca, et. al. (2003), 
in a comprehensive study which focuses its attention on the behavior of the power law index 
and QPO frequency in several BH candidate sources. We will use the data of Vignarca, our 
studies and that of other as well in our analysis (see §3.2 below). 

Observation of low term tracking of the QPO frequency with the variations of the 
photon index has been observed in at least two BH sources. Tracking is seen in the data 
of Rossi (2003) and Homan (2001) in observations of XTE J1560-500 taken over a 30-day 
period. In this period the source apparently transits from the low hard state (r ~ 1.5) to 
the intermediate or soft state where the index apparently remain close to the constant at 
a value ( T ~ 2.5 ). In this transition the QPO frequency varies from about 1 to 10 Hz. 
Similar tracking is observed in XTE J1550-564 by Sobzcak et al. (1999), (2000). However, 
the index and QPO frequency does not track with the total X-ray flux (Vignarca et al. 2003, 
Figure 7). 
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3.2. Comparison of the inferred index-frequency correlation with the data 

3.2.1. Low/ Hard- High/ Soft State Transitions 
GRS1915+105 

In Figure 7 we present the results of fit to the data of Vignarca, et. al (2003) for the 
plateau observations of the microquasar GRS 1915+105 using the TL model. The best fit 
parameters are m = 12, 5 = 1.25 and A = 1. We note that the value m = 12 is consistent 
with the BH mass evaluation m = 13.3 ± 4 obtained by Greiner et al. (2001) using the IR 
spectra and by. Shrader and Titarchuk (2003) (see also Borozdin et al. 1999) using the X-ray 
spectra. The observable index saturations at low and high values of v L are nicely reproduced 
by the model. The low and high frequency plateau regions of the data and fit are signatures 
of the two spectral phases which are explained by two different regimes of Comptonization ( 
upscattering) : (1) bulk flow Comptonization in the soft state ( saturation at T ~ 2.8) and 
(2) thermal Comptonization in the hard state ( the photon index tends to level at T ~ 1.5). 

A comparison of the index-frequency correlation of plateau observations with Vignarca’s 
(3 and v data [Belloni (1999) classification scheme] and the model fits are plotted in Figure 8. 
The only difference in the data is different saturation values at high v L that can be readily 
explained by a change of plasma temperature in the converging inflow kT C i. For upper value 
of T = 2.81 the temperature kT ci = 5 keV as for T = 2.7 the temperature kTd — 9 keV. The 
best-fit parameters for the latter data are the same as for the plateau observations except 
that in the model function r(r 0 ) (see Eq. 17) the coefficents bo, by and r* are replaced by 
bo — 2.7, bi = 22.4 and r* = 0.8 for r 0 > 2.47. 

In Figure 9 we present a plot of power-law index versus QPO centroid frequency for 
the observations of class a and u of GRS 1915+105 from Vignarca et al. (2003) (black 
points=obs. 18). Blue points correspond to the values for observations by Fiorito et al 
(2003) for OBSID 20402 and OBSID 10258. Magenta points correspond to positions where 
the OBSID 20402 data points are shifted to one half the measured QPO low frequency 
values. This suggests that the latter data points are related to the second harmonics of the 
ui frequencies. It is well known that data for NS QPOs is a mixture of the first and the 
second harmonics of horizontal branch oscillation (HBO) frequencies (see e.g. van der Klis 
2000). Red points correspond to positions where should be points with half of frequencies 
for Vignarca’ s obs. 18. A theoretical curve (blue solid line) obtained using r(r 0 ) (see Eq. 
17) (for coefficents b 0 = 2.7, bi = 22.4 and r* = 0.8 for r 0 > 2.47) for m = 12 and r 0 = 7 15 ) 
fits the Fiorito’s OBSID 10258 data points. 

It is worth noting that relation between 7 and r 0 is slightly different from that we 
obtained for the best-fits presented in Fig. 8. Namely, the index of 7, increases from 1.25 to 
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1.5. It means that for the same value of 7 the optical depth r 0 is slightly higher than that 
in state presented in Fig. 8. This change can be explained by the slightly higher plasma 
accummulation in the transition layer. 

It is interesting that the suggested half-frequency points (magenta and red ones) are 
located in a narrow corridor between these inferred curves (presented in Fig. 8 and Fig. 9 
respectively). 

XTE J1 550-564 

Figure 10 shows a plot of power law index versus QPO frequency for XTE J1550-564 
from Vignarca et al. The data again shows a flattening out or saturation of the index at 
a value T ~ 2.8. It is clearly seen from Fig. 8 of Vignarca et al. that most of the points 
presented there can be fitted by the smooth line but a few points on the right that we exclude 
for a consideration for a while. The shift transformation v' L = 12^/10 of the GRS 1915+105 
curve, T(u) (see Fig. 7) into T(v') produces a fit for the XTE 1550-564 data. This means 
that if two sources in the same accretion regime [i.e. when 70(7) is the same] their relative 
BH masses can measured by the simple shift transformation iS L = miUL/m 2 . In our case we 
use mi = 12 for GRS 1915+105 and m = 10 for XTE 1550-564. These values are very close 
to the values that have been obtained by Shrader k Titarchuk (2003) for these sources using 
X-ray spectroscopic methods (Shrader k Titarchuk 1999). Thus the TL model provides an 
independent way to estimate the mass of one or more BH’s in the same accretion state. 

Comparison of all data points of the correlation of photon index vs QPO low frequency 
for XTE 1550-564 (black points) with the inferred correlation is shown in Figure 11. The 
saturation index value in the theoretical curve is related specific value of the plasma temper- 
ature of the converging inflow kTd, which changes from 5 keV to 20 kev from the top curve 
to the bottom respectively. Slightly different cooling regimes of the site of the converging 
bulk inflow can explain the different saturation levels of the index observed in XTE 1550-564. 
The temperature of the converging inflow is determined by the disk illumination geometry 
and by the intrinsic heating of the flow by pairs produced very close to horizon (see Laurent 
k Titarchuk 2004). All these effects may be manifested by the observed index-frequency 
correlation at high values of v L . 

4U 1630-47 

Figure 12 shows a plot of index versus QPO frequency for 4U 1630-47 [Tomsick k Kaaret 
(2000), Kalemci 2002). The data again show a similar behavior to that of XTE J1550-564 
with an index saturation at T = 2.45 and F = 2.7. One can sees two phases with the index 
level at T = 1.7+0.1 and F = 2.6+0.1 and the transition between them that is much sharper 
than that observed in the other sources . The best-fit curve ) (solid line) is obtained 
using r 0 = I.67 036 and F(r) function (see formula 17) for the best-fit BH mass m = 16 
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which is close to the estimate made by Borozdin et al. (1999) using the X-ray spectroscopic 
method. The dash curve is for saturation index value T = 2.45 (related to kT C i = 20 keV). 
In this case b 0 = 2.45, bi = 53.7 and r* = 0.55 for r 0 > 2.2). Our explanation for the 
observed anomalously sharp transition is as follows. A sharp phase transition may occur 
in a source where a thick disk is formed which penetrates the coronal region as the mass 
accretion rate increases (the underlying disk is probably evaporated by corona at small mass 
accretion rates). The cool thick disk then provides an additional source for coronal cooling 
in addition to the the external disk illumination of the corona. The corona gets cold faster 
and the bulk motion Comptonization becomes a dominant process in the upscattering of soft 
disk photons for that T = 2.5 — 2.9. 


3.2.2. Very High/Soft or Extended Power Law States 

GRS1915+105 

As mentioned above observations in the “x ” states of GRS1915 by Vignarca et al. 
show a monotonically increasing linear correlation of frequency with index over a decade 
increase in QPO frequency. In this respect these results are similar to the observations 
of Kalemci (2002). What is particular interesting about Vignarca’s GRS1915 data is that 
it completely complements the data presented by Kalemci, and shows that GRS1915+105 
is quite similar to many other BH candidate sources with respect to the range of indices 
and QPO frequencies observed . FMS have also studied the behavior of the power index 
with QPO frequency in other plateau “x” states of the microquasar GRS1915+105. Like, 
Vignarca, et. al., they have employed, a simple power law and thermal disk component 
as the primary constituents of the spectral model used to fit the spectra. FMS show the 
correlations observed in the three “x” states, as mentioned above. The difference between 
the FMS results and those obtained by Vignarca et al. is that the former’s data extends 
the range of indices and QPO frequencies to higher values observed in the data analyzed by 
Vignarca et. al. In the substates observed by FMS, the positive correlation of frequency and 
index for some OBSID’s appears to be unbounded, i.e. there is no indication of saturation of 
the index value. Also the correlations observed by FMS for each of the OBSID’s is different, 
i.e. the data' does not fall on one line as is apparently the case for the Vignarca data, but 
falls on separate lines two of which appear to have the same slope. 

It is interesting to compare the extended power law index (T > 3) results mentioned 
above with those of Rau and Greiner (2003), hereafter RG. who also observe a strong position 
correlation between photon index and power law normalization observed in “x” states of GRS 
1915+105. They have found that their data fall on two distinct positive correlation curves 
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as shown in their Figure 5. Their data (like FMS) shows an extended region of index values 
well beyond the saturation values observed in the soft state , i.e. T « 2.8 up to values as high 
as T = 3.5 . This very high/soft (thermal) state is also observed in the data of Vignarca, 
et. al. for the source XTE J1748-288, where all the observed indices range from F = 2.7 — 3 
and increase monotonically with QPO low frequency. This is the same type of behavior 
observed in the extended power law data of GRS1915+105 mentioned above. 

We offer the following possible explanation for this data. The very high indices observed 
are due to a softening of the entire spectrum due to downscattering of escaping upscattered 
photons emerging from the central Comptonizing region by an extended cloud or outflowing 
wind of cold electrons well outside the TL region. This leads to an unbounded softening 
of the power law component. This type of spectral behavior has been characterized by 
Falcke (2003) as a possible jet dominated spectral state, referring to Comptonization of disk 
photons off a strong outflow of matter surrounding the BH. 

RG have also compared the relationship between radio flux and X-ray power law index. 
RG observe a positive correlation between the X-ray photon index and the radio flux, a 
negative correlation between X-ray flux and radio flux for hard X-ray photons (E > 20 keV) 
and no correlation between X-ray flux and radio flux for 1-200 keV X-rays. In terms of our 
model the observed correlation is expected for the high energy photons but not for the bulk 
photon flux in the 1-200 keV energy band which is mostly determined by the low energy 
disk photons. On the other hand the soft photons are less affected by downscattering and 
thus one does not expect correlation of the jet power (outflow rate) with 1-200 keV photon 
flux. Thus the radio data is in agreement with our model. 


3.3. Turnover States 
GRO J 1655-40 

Vignarca’s Figure 6 plots the index— frequency for two sources XTE J1550-564 and GRO 
J 1655-40 together. The behavior of the later sources is quite distinct from all others analyzed 
- it shows a reverse or negative correlation between photon index and frequency. This is the 
only case known to exhibit this type of behavior. It is also noted that the two data sets 
appears to complement one another even though they occur in separate sources. Such a 
complementarity though with different slope has been noted above between XTE J1 748-288 
and XTE Jl 550-564. It is worth noting that the saturation plateau and turnover occur at 
indices lower than 2.8-2.9. This turnover and plateau behavior can be explained in terms 
of our model either by inefficient cooling of the converging flow by the soft disk photons or 
by heating of the converging flow corona by emergent pairs generating near the BH horizon 
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(see Laurent k Titarchuk 2004 for details). Indeed, at high mass accretion the outflow pairs 
are quenched by the accreted electrons and the pairs deposit their energy in the flow. The 
increase of the converging flow plasma temperature leads to more efficient upscattering and 
thus to lower values of index (LT99). 

3.4. Correlations between High and Low Frequency QPO’s 

Strong correlations between high and low frequency QPO’s in cases where both can be 
simultaneously observed (see e.g. Swank 2001 for GRS1915+105 and Homan 2001 for XTE 
J1550-564). This type of correlations indicates the presence of nonlinear resonances which 
are expected in the inner most region near the BH. It is evidence that the region of origin 
of the QPOs is indeed very close to the BH as the TL model predicts. 

In Figure 2 we present data points showing low-high frequency dependence from BPK 
(see also PBK), Mauche (2002) and Woudt k Warner (2002) along with a fitted line u [ow = 
0.081i 'high- The plot combines data points for BH, NS and white dwarf (WD) systems. 
A reader can find the details of these observations and the QPO data analysis in BPK, 
PBK and Mauche (2002). Each data point is strickly observational, representing two QPO 
frequencies that are detected at the same time in some particular source. TW02 interpreted 
this observational correlation and identification of v hig h, vi ow with v K and vma respectively, 
lead them to determination H/Rout = 1.5 x 10 -2 and /3 = 10 (where H is the half- width of 
the disk and ft is a ratio of the gas pressure to the magnetic pressure). 


4. Summary and Conclusions 

We have developed a model of accretion onto a Black Hole which greatly simplifies and 
reclassifies the plethora of “states” observational assigned to categorize the X-ray observa- 
tions of variable BHC into two generic phases (states): 

I) soft state where we see the effect of the BH as a “drain” . Bulk inflow upscattering 
and its quenching effect on the Comptonization of disk photons dominates the behavior of 
the BH spectrum. The power law spectrum is steep in this situation. The observed high 
energy photons are emitted from a compact region, where soft energy photons of the disk are 
upscattered by bulk matter inflow forming the steep power law with photon index around 
2.8, and low QPO frequencies (above 10 Hz) and high QPO frequencies (of the order of 100 
Hz) are observed. In terms of the relativistic particle acceleration in the region of the corona 
and its effect on upscattering the soft disk photons, the soft state spectrum is a result of the 
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first order Fermi acceleration with respect to V/c. 

We would like to emphasize that bulk inflow is present in BH when the high mass 
accretion is high but not in NS, where the presence of the firm surface leads to the high 
radiation pressure which eventually stops the accretion. The bulk inflow and all its spectral 
features are absent in NSs, in particular, the saturation of the index T ~ 2.8 with respect 
to QPO frequency, which is directly related to the optical depth and mass accretion rate, 
observed in the soft spectral state of only BHs and is therefore a particular signature of a 
BH. 

II) hard state, which is comparatively starved for accretion. The hard phase (state) is 
related to an extended thermal Compton scattering cloud (cavity) characterized by a photon 
index T ~ 1.5 and the presence of low QPO frequencies (below 1 Hz). The low-hard spectrum 
is a result of the Fermi acceleration of the second order with respect to (V/c) 2 . The effect 
of the first order on V/c is smeared out by the quasi-symmetry of the particular dynamic, 
predominantly thermal motion of the Compton cloud plasma. 

The Transition Layer Componization Model: i. specifies the precise behavior of the 
power law index for each state; ii. identifies the origin and nature of low frequency and 
high frequency QPOs; iii. predicts and explains the observed relationship between them and 
their dependence mass accretion rate and spectral photon index for each of the phases; iv. 
predicts and explains the observed robust correlation of low frequency QPOs with spectral 
index observed in BHCs; v. predicts and explains the observed correlation between low and 
high frequency QPOs in BHs, neutron stars and white dwarfs. We interpret the correlation 
between low frequency v iow and power-law photon index T in quasi-periodic oscillations 
(QPO) extensively investigated by Vignarca et al. (2003) and Kalemci (2002) for a variety 
of BH sources and states and also by Homan et al. (2001) and FMS. The observed correlation 
strongly constrains theoretical models and provides clues to understanding the nature of the 
QPO phenomena particularly in BHs and in compact objects in general. 

In particular: (1) We find that the observed low frequency - the index correlation is 
a natural consequence of an adjustment of the Keplerian disk flow to the innermost sub- 
Keplerian boundary conditions near the central object. This ultimately leads to the for- 
mation of the sub-Keplerian transition layer (TL) between the adjustment radius and the 
innermost boundary (the horizon for BH). 

(2) In the framework of the TL model v high is related to the gravitational frequency at 
the outer (adjustment) radius v g re u K and vi ow is related to the magnetoacoustic oscillation 
frequency I’m a- Using a relation between Vma and the mass accretion rate and the photon 
index F and the mass accretion rate we infer a correlation between v MA and the spectral 
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index I\ 

(3) Identification of vi ow with v MA , allow us to make a comparison of the theoreti- 
cally predicted correlation with the observed correlation. For this identification we use the 
one temperature plasma assumption. We apply the plasma temperature obtained from the 
Comptonization spectra (electron temperatures) for calculations of magnetoacoustic frequen- 
cies which strongly depends on the proton temperature. The one temperature assumption 
is really consistent with the data. 

(4) We present strong arguments that in BHs the two particular distinct phases occur 
in which one of them, the steep power-law phase is a real BH signature. 

(5) We found that a hard phase (state) related to an extended Compton cloud (cavity) 
characterized by the photon index around 1.5 and the low QPO frequencies below 1 Hz. 
This is the regime where thermal Comptonization dominates the upscattering of soft disk 
photons and the spectral shape (index) is almost independent of mass accretion rate. 

(6) We find that the soft phase (state) is related to the very compact region where soft 
energy photons of the disk are up scattered forming the steep power law spectrum with 
photon index saturating around 2.8. This is the regime where bulk flow Comptonization 
dominates and the effect of an increase in the mass accretion is offset by the effect of photon 
trapping in the converging flow into the BH. 

(7) A extended power law regime where strong outflows obscure the bulk in flow region 
can also be present when the cold outflow from winds downscatters the emerging high energy 
photons and soften the observed spectrum. A regime of high accretion rate can also exist 
in this state where pair production heating adds to and can dominate the heating in the 
cavity causing a transition back to thermal Comptonization phase. In this phase a reverse 
correlation between low QPO frequency can result and is our explanation of the single case 
where this phenomena is observed. 

(8) In the context of distinction between NS and BH sources it is worth pointing to 
a certain situation that can occur in low hard states. The source should be a NS if the 
soft blackbody like component with the color temperature of order of keV is present in the 
spectrum (Torrejon et al. 2004). In contrary, such a high temperature is necessary related 
to the high disk luminosity in a BH case. 

(9) We offer a new method of BH mass estimate using the index-frequency correlations. 
Namely, if the theoretical curve of the index-frequency dependence T(ul) related to the 
parameter mi fits the data for a given source then the simple slide of the frequency axis 
v' L = (mi/m2)u L with respect to v L may allow us to obtain the mass m 2 by fit of T{v' L ) to 
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the observable correlation for another source. 

We acknowledge Tomaso Belloni for kindly supplying us data shown in Figures 7-12. 
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Fig. 1— Transition layer (corona) concept. This picture renders the gravitional energy 
release in the disk and corona and the coronal illumination by the disk soft flux in compact 
object (NS and BH). See details in text 
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High Frequency (Hz) 


Fig. 2. — Correlation between frequencies of QPO and noise components white dwarfs [dia- 
monds for SS Cyg (Mauche 2002); squares for VW Hyi (Woudt & Warner 2002)] , neutron 
star (open circles) and, black hole candidate (filled circles) sources. Neutron star and black 
hole data are from Belloni, Psaltis & van der Klis (2002). The dashed line represents the 
best-fit of the observed correlation (see TBW01 and TW02 where this correlation is predicted 
and explained using the TL model). This plot also appears in Mauche (2002). 
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Fig. 3. — Plot of QPO low frequency vs. 7 parameter (oc M). 
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Fig. 5. — Energy spectral index vs the TL optical depth. (Thermal Comptonization case) 
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Optical depth r 0 

Fig. 6. — Photon spectral index vs the TL optical depth. (General case) 
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QPO low frequency (Hz) 


Fig. 7. — Plot of power-law index versus QPO centroid frequency for the plateau observations 
of GRS 1915+105 from Vignarca et al. 2003 along with a fit using the TL model with m = 12 
and r 0 = 7 1-25 . 
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Fig. 8. — Plot of power-law index versus QPO centroid frequency for the observations of 
class ft and v (black points) and a and v (red points =obs. 15,16) of GRS 1915+105 from 
Vignarca et al. 2003 along with a fit using the TL model with m = 12 and r 0 = 7 1 ’ 25 . Values 
for the plauteau observation observations (see previous Figure) are plotted for comparison 
(blue points). 
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Fig. 9 . — Plot of power-law index versus QPO centroid frequency for the observations of 
class a and v of GRS 1915+105 from Vignarca et al. (2003) (black points=obs. 18). Blue 
points correspond to the values for observations by Fiorito et al (2003) (see text). Magenta 
points correspond to positions where should be points with half of those frequencies. Red 
points correspond to positions where should be points with half of frequencies for obs. 18 
(red points). A curve (blue solid line) is for m = 12 and r 0 = 7 1 ' 5 . 
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Fig. 10. — Comparison of the observed (Vignarca et al. 2003) and theoretical correlations 
(solid lines) of photon index vs QPO low frequency between GRS 1915+105 and XTE J1550- 
564. Black points and line for XTE J1550-564 and red points and line for GRS 1915+105. 
The XTE J1550-564 curve is produced by sliding the GRS 1915+105 curve along the fre- 
quency axis with factor 12/10 (see text for details). 
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Fig. 11. — Comparison of the observed correlation of photon index vs QPO low frequency 
for XTE 1550-564 (black points) with the inferred correlation The saturation index value in 
the theoretical curve is related to specific value of the plasma temperature of the converging 
inflow kT C i , which changes from 5 keV to 20 kev from the top curve to the bottom respectively. 
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Fig. 12. — Comparison of the observed correlation of photon index vs QPO low frequency 
with the inferred correlation. Data points are for 4U 1630-47. (see text for details). 



